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ABSTRACT 

The Perseus galaxy cluster was observed by the MAGIC Chercnkov telescope for a total effective 
time of 24.4 hr during 2008 November and December. The resulting upper limits on the 7-ray 
emission above 100 GeV are in the range of 4.6 to 7.5 x 10~ 12 cm~ 2 s _1 for spectral indices from —1.5 
to —2.5, thereby constraining the emission produced by cosmic rays, dark matter annihilations, and 
the central radio galaxy NGC 1275. Results are compatible with cosmological cluster simulations for 
the cosmic-ray-induced 7-ray emission, constraining the average cosmic ray-to-thermal pressure to 
< 4% for the cluster core region (< 8% for the entire cluster). Using simplified assumptions adopted 
in earlier work (a power-law spectrum with an index of —2.1, constant cosmic ray-to-thermal pressure 
for the peripheral cluster regions while accounting for the adiabatic contraction during the cooling flow 
formation), we would limit the ratio of cosmic ray-to-thermal energy to Ecr/E^ < 3%. Improving 
the sensitivity of this observation by a factor of about 7 will enable us to scrutinize the hadronic 
model for the Perseus radio mini-halo: a non-detection of 7-ray emission at this level implies cosmic 
ray fluxes that are too small to produce enough electrons through hadronic interactions with the 
ambient gas protons to explain the observed synchrotron emission. The upper limit also translates 
into a level of 7-ray emission from possible annihilations of the cluster dark matter (the dominant 
mass component) that is consistent with boost factors of ~ 10 4 for the typically expected dark matter 
annihilation-induced emission. Finally, the upper limits obtained for the 7-ray emission of the central 
radio galaxy NGC 1275 are consistent with the recent detection by the Fermi-LAT satellite. Due to 
the extremely large Doppler factors required for the jet, a one-zone synchrotron sclf-Compton model 
is implausible in this case. We reproduce the observed spectral energy density by using the structured 
jet (spine-layer) model which has previously been adopted to explain the high-energy emission of radio 
galaxies. 

Subject headings: galaxies: clusters: individual (Perseus) gamma rays: general 
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1. INTRODUCTION 

Clusters of galaxies provide us with the opportu- 
nity to study an "ecosystem", a volume that is a 
high-density microcosm of the rest of the Universe. 
Clusters of galaxies are the largest and most massive 
gravitationally bound systems in the Universe, with 
radii of few Mpc and total masses M ~ (10 14 — 
10 15 )M Q , of which galaxies, gas, and dark matter 
(DM) contribute r oughly for 5%, 15% and 80%, re- 
spec t ively (see, e.g.. lSarazinlll988HKochanek et al.ll2003t 
iVoitl 120051 for a general overview). While no clus- 
ter has been firm ly de t ected as a T '- ray s ource so 
far (IReimer et al.l [2001 IPerkins et all [20061: [Perkins! 



20081 lAharonian et all" 2009a. b: ]gomainko et all 120091: 



Galante et al.l 120091: iKiuchi et all l2009t lAcciari et all 
2009), they are expected to be significant 7-ray emit- 



ters on the following general grounds. (1) Clusters are 
actively evolving objects and being assembled today, in 
the latest and most energetic phase of hierarchical struc- 
ture formation. (2) Clusters serve as cosmic energy 
reservoirs for powerful sources such as radio galaxies 
and supernova-driven galactic winds. (3) Finally, clus- 
ters contain large amounts of gas with embedded mag- 
netic fields, often showing direct evidence for shocks and 
turbulence as well as relativistic particles. For recent 
reviews regarding non-thermal p rocesses in cluste rs as 
well as numerical simulations, see iBlasi et al.l (H)07) and 
iDolag et all pOOl . 

In the cosmological hierarchic clustering model, 
large-scale structures grow hierarchically through 
merging and accretion of smaller systems into larger 
ones, and clusters are th e latest and most mas- 
sive objects to form (e.g., [Peebles 1993). Recently, 
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high resolution X-ray observations by Chandra and 
XMM-Newton orbiting teles copes provided confir- 
matio n of this picture (e.g., iRosati et al.l 120021 : IVoitl 
120051) . During the course of cluster assembly, en- 
ergies of the order of the final gas binding energy 
E b — 3 x (10 61 - 10 63 ) erg should be dissipated through 
merger and accretion shocks (collectively called "struc- 
ture formation shocks") as well as turbulence. The 
energy is expected to be dissipated on a dynamical 
timescale of Td yn ~ 1 Gyr. Hence the corresponding 
rates of energy release are L ~ (10 45 — 10 47 ) erg s" 1 , 
so even a small fraction of this energy channeled 
into non-thermal particles can be of major observ- 
able consequence. Shocks and turbulence are also 
likely to accelerate non- thermal elec t rons and protons 
to hig h energies (e.g ., Uaffd 1977t ISch lickciscr et all 
19871; iBrunetti et all 120011: iMiniati et all I2001bll. " 



Ohno et all l2002t iMiniatil I2002L 120031: iSarazinl 1200 
Brunetti et all 120041: Unoue et all 12003: IBrunetti et all 
20071: IBrunetti fe Lazarianl 120071: IPfrommer et all 120071 



2008; IPfromm er 2008; Fal ceta-Goncalves et alll2010D ~ 



Clusters are also home to different types of energetic 
outflows, and the intra-cluster medium (ICM) can func- 
tion as an efficient energy reservoir. Most clusters are 
seen to harbor radio galaxies around their central re- 
gions, whose large, powerful jets of relativ istic plasma 
are interacting vigorously with the ICM (|Heinz et all 
[19981: iForman et al][20(jl iFabian et al1l2006D . A crude 
estimate of the total energy output by a single pow- 



erful radio galaxy is Erg ~ (10' 



do 



10 



<)2 i 



erg, 



tak- 



ing reasonable values for the kinetic luminosity Lrq ~ 
(10 45 — 10 46 ) erg s -1 and effective duration of a ctivity 
t RG ~ (10 7 - 10 8 ) yr (jMcNamara fc Nulsenl [2001 . The 
integrated output from the whole c luster radio galaxy 
popu l ation should be even greater ( En filin et al.l Il997l 
119981: Unoue fc Sasaki 120011) . Although rarely seen in 
present-day clusters, another source which should have 
been active in the past are galactic winds, i.e. out- 
flows driven by th e joint action of numerous supernovae 
(jVolk et al.l 119961 ). Taking the observed mass of Fe in 
the ICM to be M FeJC M ~ 3 x (10 9 - 1O 1O )M , the 
energy and Fe mass ejected by each supernovae to be re- 
spectively Esn ~ 10 51 erg and M FKi sn ~ O.lMfp, an d 
an outflow efficiency £gw ~ 0.1 (jVeilleux et al.l 12005). 
we estimate the total galactic wind energy output to be 
E G w ~ ^GwEsN/M Fe . SN M Fe j C M ~ 3 x (10 60 - 10 61 ) 
erg. In any case, along with dumping energy, these 
sources can inject substantial quantities of non-thermal 
particles into the ICM, or could have done so in the past. 

Faraday rotation measurements provide a powerful 
tool to pro be the strength of the intra-cluster mag- 
netic fie]dsjffim|etan Q99^ and even their distribu- 
tion ([Clarke et alll200*Tf ). resulting in the ICM now being 
known to be perm eated by magnetic field s with strengths 
B ~ (1 - 10) fj,G (jCarilli fc Tavlorl 127)021 : IVogt fc Enfilul 
I2005T ) . which allow for particle acceleration in shocks up 
to 7-ray emitting energies. Observations of radio ha- 
los and radio relics have already established that syn- 
chrotron emitting electrons with energies reac hing ~ 
10 GeV are present in at least some clusters (jFerettil 
120031 : iFerrari et al.l [20081 ) . although their precise origin 
is still unclear. Similar populations of electrons but with 
harder spectra may produce 7-rays efficiently via inverse 
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Compton (I C) up-scattering of the cosmic microwave 
background (lLoeb fe Wa xman 2000; Totani fc Kitavamal 
l2000riMiniatill2002ll2003tlPetrosian et al.ll2008D . Obser- 
vations in the hard X-ray regime may suggest the pres- 
ence of a non-thermal component due to the IC scatter- 
ing of cosmic microwave p hotons by relativistic electrons 
(see iRephaeli et al.1 120081 for a recent review). However, 
lAiello et al.l (|2009j) found no evidence of a hard tail above 
the thermal emission in a Swift /Burst Alert Telescope 
(BAT) sample of clusters. The ICM gas should also pro- 
vide ample target matter for inelastic collisions leading to 



pion-decay 7-rays (IVolk et alj Il996t lEnfilin et al 
Pfrommcr & Enfilml I2003L l2004at IPfrommer et al 



1997 



2008 



Pfrommcr 2008]) as well as secondary electron injection 
(iDennisonl 119801: iVestrandl Il982t IBlasi fc Colafrancescol 
1999t iDolagfc Enfflinl l2000t IPfrommer fc EnfilinT T2004a: 
Fuiita et al.ll2007t lPfrommerll2008h . The magnetic fields 



play another crucial role by confining non-thermal pro- 
tons within the cluster volume for longer than a Hub- 
ble time, i.e. any protons injected int o the ICM accu- 
mulates throughout th e cluster history (|Volk et al.lll996l : 
iBerezinskv et al.lll997D . 

Galaxy clusters present very large M/L ratios and con- 
siderable overdensities, which are crucial for indirect DM 
searches. Despite the fact that they are not as near as 
other po tential DM candidates, as the dwarf spheroidal 
galaxies (|Albert et al.ll2008dlAliu et al.l2009al) . the large 
DM masses of clusters could make them ideal laborato- 
ries also for the sear c h of a DM annihila tion 7-ray signal 
(|Jeltema et al.ll2009l: IFinzke et al.H2009[) . 

In this paper we report the results of the Perseus clus- 
ter observation performed by the Major Atmospheric 
Gamma Imaging Cherenkov (MAGIC) telescope for a to- 
tal effective time of 24.4 hr during 2008 November and 
December. In Section 2, we explain the physical moti- 
vations why we chose Perseus over other galaxy clusters 
and present its main characteristics. In Section 3, we 
briefly introduce the MAGIC telescope. We then de- 
scribe the Perseus data sample, the analysis, and the 
obtained flux upper limits. We discuss the implica- 
tions for the cosmic ray (CR) pressure and the possi- 
ble DM annihilation-induced 7-ray emission in Section 4 
and 5, respectively. In Section 6, we discuss the impli- 
cations for the jet emission model of the central radio 
galaxy NGC 1275. Finally, in Section 7, we summarize 
our conclusions. All cluster masses and luminosities are 
scaled to the currently favored value of Hubble's constant 
H = 70 km s- 1 Mpc" 1 . 

2. TARGET SELECTION AND PRELIMINARIES 

The Perseus cluster, also called A426, is at a distance 
of 7 7.7 Mpc (z = 0.0 18). It is the brightest X-ray clus- 
ter (jEdge et al.l Il992f ) and hosts a massive cooling flow 
with high central gas densities of 0.05 cm -3 (see Tabled]). 
Perseus furthermore hosts a luminous radio mini-halo - 
diffuse synchrotron emission that fills a large fraction of 
the cluster core region - and sh ows a source extension of 
~ 200 kpc (|Pedlar et al.lll990l) . This radio mini-halo is 
well modeled by the hadronic scenario where the radio 
emitting electrons are produced in hadronic CR proton 
interactions with ambient gas protons requiring only a 
very modest fraction of a few percent CR pressure re l- 
ative to thermal pressure (|Pfrommer fc Enfilinl l2004aD . 
In particular, the similarity of the thermal X-ray emis- 



sion to that of the radio mini-halo comes about natu- 
rally as both processes scale with the number density 
squared. An alternati ve model for the radio emission 
has been proposed by iGitti et ail (|2002| ) which explains 
the radio mini-halo by re-acceleration of relativistic elec- 
trons through second-order interactions with magneto- 
hydrodynamic (MHD) turbulence. However, it remains 
to be shown whether the necessary turbulent energy den- 
sity can be provided throughout the entire cooling flow 
region of Perseus. These conditions provide high target 
densities for hadronic CRp-p interactions and enhance 
the resulting 7-ray flux. 

The Perseus galaxy cluster was carefully chosen over 
other nearby clusters after considering the expected 7- 
ray emission from the pion-decay and DM annihilation. 
Moreover, the central radio galaxy NGC 1275 is expected 
to be a promising GeV-TeV target, and hence is another 
strong motivation to observe this cluster. In the following 
subsections, we detail our considerations. 

2.1. Cosmic-Ray-Induced Emission 

In the course of this work, we used cosmological sim- 
ulations of the formation of galaxy clusters to inform us 
about the expected spatial and spectral characteristics 
of the CR induced 7-ray emission. A clear detection of 
the IC emission from shock-accelerated CR electrons will 
be challenging for Imaging Atmospheric Cherenkov Tele- 
scopes (IACTs) due to the large angular extent of these 
accretion shocks that subtend solid angles corresponding 
to up to six virial radii. For these instruments, the spa- 
tially concentrated pion-decay 7-ray emission resulting 
from hadronic C R interactions that dominates the tota l 
7-ray luminosity (jPfrommer et al . 2008; Pfro mmerll2008[ ) 
should be more readily detectable than the emission from 
the outer region. 

To address the question of universality and predictabil- 
ity of the expected 7-ray emission, we simulated a sample 
of 14 galaxy clusters that span one and a half decades 
in mass and show a variety of dynamical states ranging 
from relaxed cool core clusters to violent merging clusters 
(details are given in Sect. I4.1j) . In order to find the most 
promising target cluster in the local Universe for detect- 
ing the pion-decay emission, we computed the scaling 
relations betwe en the 7-ray lum inosity and cluster mass 
of our sample (jPfrommer! 120081 ) and used these to nor- 
malize the CR-induced emission of all clusters in a com- 
plete sample of the X- ray brightest clusters (the ex tended 
HIFLUGCS catalog: IReiprich fc Bohrhigerll2002T) . This 
favors a high-mass, nearby galaxy cluster with a scaling 

-^foo/Aum' where M 2 oo is the virial mass 1 , Aum is the 
luminosity distance, and f3 ~ 1.32 is a weakly model- 
dependent scaling parameter that provides the rank or- 
deri ng according to t he brightness of each individual clus- 
ter (|Pfrommerll2008[ ) . As a second criterion, we required 
low zenith angle observations, i.e below 35°, that en- 
sure the lowest possible energy thresholds and the max- 
imum sensitivity for the detector. We carefully mod- 
eled the most promising targets, accounting for the mea- 
sured gas density and temperatures from thermal X-ray 
measurements while assuming a constant CR-to-thermal 

1 We define the virial mass and the virial radius Ra as 

the mass and radius of a sphere enclosing a mean density that is 
A = 200 times the critical density of the Universe. 
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TABLE l 

Properties of the Perseus Galaxy Cluster 



Aum [Mpc] ii 20 o [Mpc] M 200 [M©] i x ,o.i-2.4 [ergs" 1 ] T x [keV] L V=1A [ergs" 1 Hz" 



0.0183 



77.7 



1.9 



7.71 x 10 1 



8.31 x 10 4 



6.8 



3.38 x 10 31 



Notes. Data taken from IReiprich Hi BohriEierl J2002I) . [Pecllar et alj 1 119901) . and lChurazov et al.l (|20"03T) . 



gas ratio (|Pfrommer fc EnfiTmll2004al) . Cluster- wide ex- 
tended radio synchrotron emission that informs about 
present high-energy processes was additionally taken into 
account before we selected the Perseus cluster as our 
most promising source. Although other clusters showed 
a somewhat higher 7-ray flux in our simulations (e.g., 
Ophiuchus), the facts that Perseus is observable at low 
zenith angles and that the expected emission is more spa- 
tially concentrated make it the best suited target for this 
observation. 

2.2. Dark Matter Content 

Typically up to 80% of the total mass of a galaxy clus- 
ter is in the form of non-baryonic DM. Since the DM 
annihilation 7-ray signal is expected to be proportional 
to the integrated squared DM density along the lin e of 
sight (|Evans et alJ!2004t iBergstrom fc Hooper! I2006f >. it 
is obvious that galaxy clusters could be good candidates 
to look for DM as well. This is true despite the fact that 
they are located at much larger distances than other po- 
tential DM candidates, such as dwarf spheroidal galaxy 
satellites of the Milky Way or the Galactic Center. One 
obvious reason is the huge amount of DM hosted by clus- 
ters compared with the rest of candidates. Perseus, for 
example, is located ~1000 times farther than Milky Way 
dwarfs, but it contains roughly six orders of magnitude 
more DM than the Willman 1 dwarf galaxy, one of the 
most promising DM candidates accord ing to recent work 
(jStrigari et al.ll2007t lAliu et al.ll2009a[) . Additionally, the 
presence of substructures could be of crucial importance. 
Substructures in clusters may significantly enhance the 
DM signal over the smooth halo, while we do not expect 
this to be of special relevance for dwarf galaxies since 
their outer regions are s everely affected by tidal strip- 
ping (jPinzke et al.l[2009t Sanchez-Conde et al. 2010, in 
prep.). 

Essentially, the annihilation flux is pro portional to 
the p roduct of two parameters (see e.g., Evans et al.l 
120041 for details): a first one that captures all the par- 
ticle physics (DM particle mass, cross section, etc.), 
which we will label as /susy, and a second one, J as tro, 
that accounts for all the astrophysical considerations 
(DM distribution, telescope point-spread function (PSF), 
etc.). The particle physics factor just acts as a nor- 
malization in the expected annihilation flux, so we can 
neglect it when performing a comparative study - as 
we are doing in this section. Concerning the astro- 
physical factor, the DM distribution is commonly mod- 
eled with radial density profiles of the form p(r) = 
p s /[(r/r s ) 7 (1 + {r/r s ) a )^~^/ a ], where p s and r s rep- 
resent a characteristic densi ty and a scale radius respec- 
tively (jKravtsov et al.l ll998). These density profiles are 
well motivated by high-resolution N-body cosmological 
simulations. Here we adopt the Navarro-Frenk- White 



(jNavarro et al.l 119971 hereafter NFW) DM density pro- 
file, with (a,/3,j) = (1, 3, 1). For an NFW profile, 
90% of the DM annihilation flux comes from the region 
within r s , so that the corresponding integrated lumi- 
nosity is proportional to r^p 2 s . We can derive r s and 
p s for Perseus, assuming M 2 oo = 7.7 x 10 14 M Q (as 
given in Table Ht and a concentration of ~6 (as given 
by the lBullock et al]|2001l virial mass-concentration scal- 
ing relation). We obtain r s = 0.384 Mpc and p s — 
1.06 x 10 15 M Q Mpc -3 , which translates into a total value 
of Jastro ~ 1.4 x 10 16 GeV 2 cm -5 for the scale radius re- 
gion. In the case of Coma, although slightly (~15%) 
more massive than Perseus, the fact that it is located 
significantly farther (101 Mpc) translates into a slightly 
lower annihilation flux. Virgo, only 17 Mpc away from 
us, gives a larger DM annihilation flux, but here the large 
extension of the region from which most of the annihi- 
lation flux is expected to come compared with Perseus 
(r s ~ 1°.2 and r s ~ 0°.3, respectively) could represent an 
obstacle from the observational point of view. Source ex- 
tension is of special relevance for single-telescope IACTs, 
for which point-like sources (sources with an angular ex- 
tension smaller than or similar to the telescope PSF) are 
more readily observable. 

2.3. The NGC 1275 Radio Galaxy 

The central NGC 1275 radio galaxy is another strong 
motivation for 7-ray observations of the Perseus galaxy 
cluster. The detection at TeV ene rgies of the radio 
;alaxies M 87 (lAharonian et al.l 120061 ) and Centaurus A 
Aharonian et al. 2009c) has forced a substantial revision 
of the paradigm whereby very high energy (VHE) emis- 
sion is a characteristic property of highly relativistic jets 
closely aligned with the line of sight, establishing radio 
galaxies as a new class of VHE 7-ray emitters. Note that 
NGC 1275 has various characteristics in common with 
Centaurus A which has also been inte rpreted as a pos- 
sible source of ultra-high-energy CRs (|Hardcastle et al.l 
120091 ). 

The NGC 1275 radio galaxy is the brightest radio 
source in the northern sky. Its jet inclination an- 
gle seems to increase from 10° — 20° at m illiarcsecond 
scale s up to 40° — 60° at arcsecond scales (|Dunn et al.1 
12006ft. Note that NGC 1275 was classified as a blazar 
by"! ngel fc Stockman! (|1980f ) because of its optical po- 
larization, and it has be en seen to vary in the optical 
on time scales of a day (|Geller et al.l 11979ft . All these 
elements are promising from the point of view of the 
TeV detectability, since they suggest that the emission 
region is located at the base of the jet. In these condi- 
tions, in the scenario based o n the structured jet model 
(jTavecchio fc Ghiselhnl I2008D . we expected VHE emis- 
sion from the layer of the jet at a level detectable by 
MAGIC. 
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3. MAGIC OBSERVATION AND RESULTS 

The MAGIC telescope is located on the Canary Island 
of La Palma (2200 m asl, 28.45°N, 17.54°W). With a pri- 
mary mirror diameter of 17 m, it is currently the largest 
IACT. CRs impinging the Earth atmosphere originate 
atmospheric showers that in turn produce Cherenkov 
light. The ultra-violet Cherenkov flashes are reflected in 
the focal plane of the telescope, where a camera of 577 
photomultipliers records the resulting images. MAGIC 
reconstructs the incoming 7-ray directions with an ac- 
curacy of about 0°.l and achieves an energy resolution 
above 150 GeV o f about 20% (see lAlbert et all l2008d: 
lAliu et al.ir2009bl for details). 

3.1. Observation and Analysis 

MAGIC observed the Perseus cluster for 33.4 hr during 
2008 November and December, at zenith angles between 
12° and 32°, which guarantees the lowest energy thresh- 
old. The observation wa s performed in the false-source 
tracking (wobble) mode ([Fomin et al.lll994T) pointing al- 
ternatively to two different sky directions, each at a 0°.4 
distance from the nominal target position. 

The main background for Cherenkov telescopes is due 
to the hadronic CRs and the night sky background. Our 
standard analysi s procedure is as fol lows (for a detailed 
description, see lAlbert et al.l 12008c) : data calibration 
and extra ction of the number o f photoelectrons per pixel 
are done (| Albert et al.l l20"08aD . This is followed by an 
image cleaning procedure using the amplitude and tim- 
ing information of the calibrated signals. Particularly, 
the arrival times in pixels containing > 6 photoelectrons 
(core pixels) are required to be within a time window 
of 4.5 ns and for pixels containing > 3 photoelectrons 
(boundary pixels) within a time window of 1.5 ns from 
a neighboring core pixel. For the surviving pixels of 
each event, the shower parameters are reco nstructed us - 
ing the Hillas parameterization algorithm (Hillas 1985). 
Hadronic background suppression is achieved using a 
multivariate method cal led Random Forest (jBreimanl 
[200lHAlbert et al.|[2008bl) . which uses the Hillas param- 
eters to define an estimator called hadronness (it runs 
from for gammas to 1 for hadrons) by comparison 
with Monte Carlo (MC) 7-ray simulations. Moreover, 
the Random Forest method is used for the energy esti- 
mation of a reconstructed shower. The gamma/hadron 
(g/h) separation in the analysis was optimized on a sam- 
ple of well-understood Crab Nebula data, which is com- 
monly accepted as a standard reference source for VHE 
astronomy. 

Part of the data has been rejected mainly due to the 
bad weather conditions during some observation days. 
The total data rejected amount to ~ 27%, resulting in 
24.4 hr effective observation time of very high data qual- 
ity. Independent cross-checks were performed on the 
data giving compatible results. 

3.2. Results 

Given the good data quality and the low zenith angles 
of observation, the analysis energy threshold results to 
be 80 GeV. Beyond this threshold, no significant excess 
of 7-rays above the background was detected in 24.4 hr 
of observation. In Figure [TJ the a-plot for energies above 
250 GeV, where the best integral sensitivity is obtained 



from a Crab Nebula data sample, is reported. The a- 
parameter is denned as the angular distance between 
the shower image main axis and the line connecting the 
observed source position in the camera and the image 
barycenter. Background events are isotropic in nature, 
and thus produce randomly oriented shower images. This 
results in a more or less smooth event distribution in the 
a-plot. The 7-ray events due to the source, on the other 
hand, are predominantly aligned to the observed posi- 
tion in the camera. For a detected source, this results in 
a significant excess of events at small a. A fiducial re- 
gion a < 6° and a hadronness cut of 0.05 are chosen by 
optimizing the analysis on a Crab Nebula data sample. 
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Fig. 1. — Perseus a-plot as seen by MAGIC in 24.4 hr above 
250 GeV using a hadronness cut of < 0.05. The blue crosses rep- 
resent the signal, and the red shaded region is the background. 
The vertical black dotted line represents the fiducial region a < 6° 
where the signal is expected. Only events above 250 GeV are dis- 
played since the best integral sensitivity, around 1.6% of Crab, is 
obtained from a Crab Nebula data sample in this energy range. 

In Figure [21 the significance map for events above 
150 GeV in the observed sky region is shown. The source- 
independent DISP method has been used. This implies 
the rise of th e energy threshold from 80 GeV t o around 
150 GeV (see iDomi ngo-San tamaria et al.ll2005l for a de- 
tailed description). The significance distribution in the 
map is consistent with background fluctuations. In Fig- 
ure [21 Xj^Y_contours from the XMM-Newton observa- 
tions (jChurazov et al.l l2003) are also shown. 

The significan ce was calculated according to Eqn. 17 of 
ILi fc Mai (Il98l and uppe r limit estimation is performed 
using the Rolke method (|Rolke et al.l [2005). The upper 
limits in number of excess events are calculated with a 
confidence level of 95%. For the upper limit calculation, 
a systematic uncertainty of 30% in the energy estima- 
tion and effective area calculation is taken into account. 
Our systematic error budget is obtained by adding up 
the individual contributions in quadrature. The different 
sources of systematic uncertainties are mainly related to 
the differences between the real experimental co nditions 
and the simulated ones ( see lAlbert et al.|[2008d for a de- 
tailed discussion on the systematic errors). The photon 
flux upper limit is finally rec onstructed for a g eneral 7- 
ray spectrum as described in lAliu et al.l (|2009al) . 

In sections [4] and [U we will discuss the implications of 
this observation for the CR and DM annihilation-induced 
7-ray flux, respectively. Using the true density profile a s 
obtained by X-ray measurements (jChurazov et alfc OOS) , 
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TABLE 2 

Integral Flux Upper Limits Above 100 GeV 



3.26 3.28 3.3 3.32 3.34 3.36 3.38 3.4 

RA [h] 

Fig. 2. — Significance map for events above 150 GeV in the ob- 
served Perseus cluster sky region. The significance distribution 
is consistent with background fluctuations. Blac k contours from 
XMM-Newton observations in the X-ray band (Chura zov et all 
2003) are also shown. The angular extent of the outermost con- 
tours is approximately 0°.45, which corresponds to ~ 610 kpc. 

we will be able to model the spatial characteristics of the 
CR-induced 7-ray signal. Our simulations indicate that 
60% of the total 7-ray flux are contained within a circle 
of radius rn.6 = 0°.15 (this angular scale corresponds to 
a physical radius of 200 kpc) . We then compare the flux 
from within this region to the upper limits. As the char- 
acteristics of the considered emission region are close to 
a point source we use point-like upper limits. The same 
conclusion is valid also for the DM annihilation signal. 
In this case, as explained in section |2~2| the 90% of the 
expected emission is coming form the scale radius region. 
For Perseus, we obtained r s ~ 0°.3 which is somewhat 
extended compared to the telescope angular resolution. 
However, the fact that the NFW profile is very steep im- 
plies that the main DM emission comes from the core of 
the source that can be considered approximately point- 
like compared to our angular resolution. 

To compute flux upper limits, we assume specific spec- 
tral indices that have been motivated by an astrophysi- 
cal scenario in mind (see the following sections). This 
"scenario-guided" approach allows us to provide the 
tightest limits on physically motivated parameters and 
underlying astrophysical models. In the next sections we 
will consider flux upper limits computed using a power- 
law 7-ray spectrum with spectral indices T of— 1.5, —2.2, 
and —2.5. In Tabled the corresponding integral flux up- 
per limits for energies above 100 GeV are listed. 

In Section |4l we will use an integral flux upper limit 
set above given energy thresholds in order to trace the 
energy range where we can better constrain the models. 
In Table [3J the obtained integral flux upper limits for 
T = —2.2 are shown. Note that we do not compute in- 
tegral upper limits above 80 GeV (as we have not shown 
a cumulative a-plot for energies above this value). This 
is because the g/h separation for events below 100 GeV 
works in a substantially different way with respect to the 
higher energy events. Therefore, we analyze separately 
the events below 100 GeV and the events of higher en- 
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Notes. Integral flux upper limits are listed for a power-law 7-ray 
spectrum with spectral index T for energies above 100 GeV. The 
corresponding upper limit for the number of excess events is 186. 

ergy, with different sets of analysis cuts. 

Finally, for completeness, in Table |4] the differential 
flux upper limits for the assumed spectral indices are 
shown in different energy intervals. Spectral energy den- 
sity (SED) upper limits can also be obtained from those 
differential flux upper limits, as done in Section [6] dis- 
cussing the observation implications for the radio galaxy 
NGC 1275. 
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Integral Flux Upper Limits for a Power-law 7-RAY 
Spectrum with Spectral Index T = —2.2 Above a Given 
Energy Threshold E th . 
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3.3. Comparison to Previous Observations 

There are fe w existing IACT o bservatio ns 
galaxy clusters (IPerkins et al.l [20061 IPerkins! 



of 



lAharonian et al.l 



Domainko et al.l 



2008 



2009; 



2009aH_ _ 

Galante et al.l 120091: iKiuchi et al.l 120091 ; lAcciari et al.l 
20091) . In section 14.31 we will compare the limits on 



the CR-to-thermal pressure obtained by other IACTs 
with those derived in this work. However, there are 
two observations of the Perseus gal axy cluster made 
by WHIPPLE (IPerkins et al.l [2006h and VERITAS 
(|Acciari et al.ll2009f) with which we can directly compare 
our upper limits. 

The WHIPP LE Collaboration o bserved the Perseus 
galaxy cluster (|Perkins et al.l [2006h for ~ 13 hr obtain- 
ing an integral upper limit above 400 GeV of 4.53 x 
10~ 12 cm~ 2 s _1 assuming a spectral index V = —2.1. 
We can compare this value with our integral upper limit 
above 400 GeV of 1.83 x 10~ 12 cm~ 2 s" 1 with V = -2.2 
(see Table [3]). Our upper limit is significantly lower 
than the WHIPPLE one; clearly, this is not a surprise 
as the MAGIC telescope belongs to a new generation 
of IACTs. More rece ntly, the VERITA S Collabora- 
tion observed Perseus (|Acciari et al.l I200"9l) for ~ 8 hr 
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Differential Flux Upper Limits 
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Notes. Differential ffux upper limits are listed in units of 10 
in energy ranges in units of GeV. 

and obtained an integral upper limit above 126 GeV of 
1.27 x 10 -11 cm~ 2 s _1 assuming T — —2.5. We can 
compare this value with our corresponding integral up- 
per limit above 100 GeV of 7.52 x 10~ 12 cm~ 2 s _1 (see 
Table Despite the fact tha t the VERITAS s ensitivity 
of about 1% of Crab Nebula (|Otte et al.ll2009f ) is better 
than the MAGIC one, our upper limi t is slightly lower 
than that found bv lAcciari et al.1 (|2009f ) as expected from 
the significant difference in observation time. 

4. COSMIC-RAY-INDUCED EMISSION 

We use the upper limits on the integrated flux (Tabled 
to put constraints on the CR-to-thermal pressure distri- 
bution and pursue three different approaches. (1) We 
perform high-resolution hydrodynamical simulations of 
cluster formation and evolution in a cosmological frame- 
work that include CR physics to predict the 7-ray emis- 
sion and to o btain limits on the CR-to-th ermal pressure. 
(2) Following IPfrommer fc Enfflinl (|2004aD . we use a sim- 
plified approach that assumes a constant CR-to-thermal 
energy density, a power-law spectrum in momentum, and 
compare the resulting CR-to-thermal pressure limits to 
those obtained by other IACT observations. (3) We use 
the observed luminosity of the radio mini-halo to place 
a lower limit on the expected 7-ray flux in the hadronic 
model of the radio mini-halo. This translates into a min- 
imum CR pressure that is crucial for disentangling the 
emission mechanism in the radio and provides a clear 
prediction for the expected 7-ray flux. 

Before doing so, we detail our cosmological simulations 
that we base our main analysis on. To this end, we in- 
vestigated the spatial and spectral properties of 7-ray 
emission in these simulations and refer the reader to the 
theory papers fo r further details ([Pfrommer et al.l 120081 : 
Pfrommcr 2008, Pinzkc & Pfrommer, in prep.). 

4.1. Cosmological Simulations 

Simulations were performed using the "concordance" 
cosmological cold DM model with a cosmological con- 
stant (ACDM) motivated by First Year cosmological 
constraints of Wilkinson Microwave Anisotropy Probe 
(WMAP). The cosmological parameters of our model are: 

n m = a dm + n h = 0.3, n h = 0.039, n A = 0.7, h = 0.7, 

n = 1, and er 8 = 0.9. Here, fi m denotes the total mat- 
ter density in units of the critical density for geometri- 
cal closure today, p a a{z = 0) = 3H^/(8ttG). fib and 
fi a denote the densities of baryons and the cosmological 
constant at the present day, respectively. The spectral 
index of the primordial power-spectrum is denoted by n, 
and us is the rms linear mass fluctuation within a sphere 
of radius 8 h~ 1 Mpc extrapolated to z = 0. 

Our simulations were carried out with an updated 
and extended version of the distributed-memory par- 



cm 2 s 1 TcV 1 for a power-law 7-ray spectrum with spectral index T 



allel TreeSPH code GADGET-2 (ISnringel et al.l 120011 
Springcl 2005). Gravitational forces were computed us- 
ing a combination of particle-mesh and tree algorithms. 
Hydrodynamic forces were computed with a variant of 
the smoothed particle hydrodynamics (SPH) algorithm 
that conserves energy and entro py where appropriate, 
i.e. o utside of shocked regions ((Springel fc Hernquistl 
l2002f K We have performed high-resolution hydrodynamic 
simulations of a sample of galaxy clusters that span over 
one and a half decades in mass and show a variety of 
dynamical states ranging from relaxed cool core clusters 
to violent merging clusters. Our simulated clusters have 
originally b een selected from a low-resolution DM-only 
simulation (Yoshida et al. 2001fl. Using the 'zo omed ini- 
tial conditions' technique ()Katz fc Wh ite 1993), the clus- 
ters have been re-simulated with higher mass and force 
resolution. In high-resolution regions, the DM particles 
had masses of md m = 1.61 X 10 9 M Q and SPH parti- 
cles m gas = 2.4 x 10 8 h^} M so each individual cluster 
is resolved by 8 x 10 4 -4 x 10 6 particles, depending on its 
final mass. The SPH densities were computed from 48 
neighbours, allowing the SPH smoothing length to drop 
at most to half of the value of the gravitational soften- 
ing length of the gas particles. This choice of the SPH 
smoothing length leads to our minimum gas resolution of 
approximately 1.1 x 10 10 hj^ M . For the initial redshift 
we chose 1 + Zj n i t = 60. T he gravitational force soft ening 
was of a spline form (e.g.. iHernquist &: Katdll989l ) with 
a Plummer-equivalent softening length that is assumed 
to have a constant comoving scale down to z = 5, and 
a constant value of 7/i 70 1 kpc in physical units at later 
epochs. 

These simulations included radiative hydrodynamics, 
star formation, supernova feedback and followed CR 
physics using a novel formulation that followed the most 
important injection and loss processes self-consistently 
while accou nting for the CR press u re in the equations 
of motion (IPfrommer et alj 120061 : lEnfilin et ail 120071: 
Uubelgas et al.ll2008f ). To obtain predictions of the GeV- 
TeV 7-ray emission from clusters, we used an updated 
version of the CR physics in our code. It is capable of 
following the spectral evolution of the CR distribution 
function by tracking multiple CR populations in each 
gaseous fluid element; each of these populations is de- 
scribed by an amplitude, a low-momentum cut-off, and a 
characteristic power-law distribution in particle momen- 
tum with a distinctive slope that is determined by the 
acceleration process at formation shocks or supernova 
remnants (Pinzke & Pfrommer 2010, in prep.). Adia- 
batic CR transport processes such as compression and 
rarefaction, and a number of physical source and sink 
terms which modify the CR pressure of each particle 
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are modeled. The most important sources considered 
are diffusive shock acceleration at cosmological structure 
formation shocks and optional injection by supernovae 
while the primary sinks are thermalization by Coulomb 
interactions and catastrophic losses by hadronic inter- 
actions. We note that the overall normalization of the 
CR distribution scales with the maximum acceleration 
efficiency at structure formation shock wav es. Following 
recen t observations at supernova re mnants dHelder et al.l 
|2009() as well as theoretical studies (|Kang fc Jonesl ll2005'). 
we adopt a realistic value of this parameter and assume 
that 50% of the dissipated energy at strong shocks is in- 
jected into CRs while this efficienc y rapidly decreases for 
weaker shocks (|Enfilin et al.l 120071 ). 

We computed the 7-ray emission signal and found that 
it obeys a universal spectrum and spatial distribution 
(Pinzke & Pfrommer 2010, in prep.). This is inherited 
from the universal concave spectrum of CRs in galaxy 
clusters that is caused by the functional form and redshift 
dependence of the Mach number distribution of structure 
formati on shocks that are resp onsible for the acceleration 
of CRs (|Pfrommer et al.ll2006f ). The CR distribution has 
a spectral index of T ~ —2.5 at GeV energies and expe- 
riences a flattening towards higher energies resulting in 
T ~ —2.2 at energies above a few TeV. Hence, the result- 
ing 7-ray spectrum from CR induced pion-decay shows 
a characteristic spectral index of V ~ —2.2 in the energy 
regime ranging from 100 GeV to TeV. The spatial dis- 
tribution of the CR number dens ity is mainly governed 
by adiabatic transport processes (jPfrommer et al.ll2007| ) 
and similarly attains an approximate universal shape rel- 
ative to that of the gas density. These findings allow us 
to reliably model the CR signal from nearby galaxy clus- 
ters using their true density profiles as obtained by X-ray 
measurements that we map onto our simulated density 
profiles. 

In addition to CR protons, we modeled relativistic elec- 
trons that have been accelerated at cosmological struc- 
ture formation shocks (primary CR electrons) and those 
that have been produced in hadronic interactions of 
CRs with ambient gas protons (secondary CR electrons) . 
Both populations of CR electrons contribute to the 7-ray 
emission through Compton up-scattering photons from 
the cosmic microwave background as well as the cumula- 
tive star light from galaxies. It turns out that the pion- 
decay emission of the cluster dominates over the con- 
tribution from b oth IC compone nts - in particular, for 
relaxed systems (|Pfrommerll2008f) . 

In our optimistic CR model (radiative physics with 
galaxies), we calculated the cluster total 7-ray flux within 
a given solid angle. In contrast, we cut the emission from 
individual galaxies and compact galactic-sized objects in 
our more conservative model (radiative physics without 
galaxies). In short, the ICM is a multiphase medium con- 
sisting of a hot phase which attained its entropy through 
structure formation shock waves dissipating gravitational 
energy associated with hierarchical clustering into ther- 
mal energy. The dense, cold phase consists of the true 
interstellar medium (ISM) within galaxies and at the 
cluster center as well as the ram-pressure-stripped ISM. 
These cold dense gas clumps dissociate incompletely in 
the ICM due to insufficient numerical resolution as well 
as so far incompletely understood physical properties of 
the cluster plasma. All of these phases contribute to the 



& 10" 



10" 





1 1 1 1 1 












radiative physics w/o gal. x 2 
















radiative physics w/o gal. 








F ymiD <5 =10^G, e s < Efc/3) 








i : 




" - „ t 


: • 



















A 

* * 




100 




1000 



E [GeV] 

Fig. 3. — Integral flux upper limits (this work, Table [3} are com- 
pared with simulated integrated spectra of the 7-ray emission from 
decaying neutral pions that result from hadronic CR interactions 
with the ambient gas in the Perseus cluster. Our conservative 
model without galaxies (solid) is contrasted to our model with 
galaxies (dashed). We scaled our conservative model with a factor 
of 2 so that it is just consistent with the upper limits obtained in 
this work (dotted). In our simulations, we assume an observation- 
ally motivated large value for the maximum CR energy injection 
efficiency at structure formation shocks and convert half of the 
dissipated energy to CRs at strong shocks. Smaller values would 
imply smaller 7-ray fluxes. Additionally shown are minimum 7-ray 
flux estimates for the hadronic model of the radio mini-halo of the 
Perseus cluster (dash-dotted with minimum flux arrows; see the 
main text for details). Note that a non-detection of 7-rays at this 
level seriously challenges the hadronic model. 

7-ray emission from a cluster. To assess the bias asso- 
ciated with this issue, we performed our analysis with 
both limiting cases bracketing the realistic case. 

In Figure [3j we compare the integral flux upper limits 
obtained in this work (see Table with the simulated 
flux that is emitted within a circle of radius ro.6 = 0°.15 
for our two models, with and without galaxies. The up- 
per limits are a factor of 2 larger than our conservative 
model and a factor of 1.5 larger than our most optimistic 
model predictions implying consistency with our cosmo- 
logical cluster simulations. We note however that our 
simulated flux represents a theoretical upper limit of the 
expected 7-ray flux from structure formation CRs; lower- 
ing the maximum acceleration efficiency would decrease 
the CR number density as well as the resulting 7-ray 
emission. 

4.2. Constraints on the Cosmic Ray Pressure 

In Figure [H we show the simulated 7-ray surface 
brightness map of a cooling flow cluster of mass simi- 
lar to Perseus. As the CR-induced 7-ray flux is a radi- 
ally declining function, so is the CR pressure. A quan- 
tity that is of great theoretical interest is the CR pres- 
sure relative to the thermal pressure Vcr = PcR./-Pth> 
as it directly assesses the CR bias of hydrostatic cluster 
masses since the CR pressure enters in the equation of 
motion. On the right-hand side of Figure [4j we show the 
profile of the CR-to-thermal pressure (volume- weighted) 
of this simulated cluster. Moving from the periphery 
towards the center, this quantity is a steadily declin- 
ing function until we approach the cooling flow region 
around the cD galaxy of this cluster (similar to NGC 
1275) where the CR pressure rises dramatically rela- 
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Fig. 4. — Left: simulated 7-ray em ission at energies E > 100 GeV from a cluster that has twice the mass as Perseus (using the simulation 
of the cooling flow cluster g51 from Pfrommcr ct al. 2008). We show the sum of pion-decay-induced 7-rays (which dominates the central 
and the total flux) and the IC emission of CR electrons accelerated at formation shocks and by hadronic CR interactions. Right: profile 
of the CR-to-thermal pressure (volume-weighted) of this cluster. We contrast a simulation where we only accelerate CRs at structure 
formation shocks of the entire cosmic history (solid) with one where we additionally account for CRs that are injected through supernova 
feedback within the star- forming regions in our simulation (dashed). 



tive to tha t of the thermal gas w hich cools on a short 
time scale (jPfrommer et al. 200f|). The volume average 
is (X CR ) = (P C R)/(Pth> = 0.02, dominated by the re- 
gion around the virial radius, while the ratio of CR-to- 
thermal energy is given by Eck/E^ — 0.032 2 . Perseus 
has a smaller mass and a corresponding temperature that 
is only half of that of our simulated cooling flow cluster. 
Noting that Xqr oc l/Ph oc 1/kT, 3 we expect these 
values to be a factor of ~ 2 larger in Perseus, yielding 
(X C r) ^ 0.04 for the entire cluster and (X C r) ^ 0.02 for 
the core region that we probe with the present observa- 
tion. 

We have to scale our conservative model prediction by 
a factor of ~ 2 to reach the upper limits (cf. Figure [3]) 
which implies that this work constrains the relative pres- 
sure contained in CRs to < 8% for the entire cluster and 
to < 4% for the cluster core region. The presence of dense 
gas clumps potentially biases the simulated 7-ray flux 
high and hence the inferred limits on Xcr low. Another 
source of bias could be unresolved point sources inside 
the cluster such as an active galactic nucleus (AGN). In 
the presented simulation of the cool core cluster g51, the 
bias due to subclumps amounts to a factor of 1.5 but it 
could be as high as 2.4 which is the mean difference be- 
tween our conservative and optimistic model across our 
scaling relations. We note however that the latter case is 
already excluded by our upper limits provided the max- 
imum shock acceleration efficiency is indeed as high as 
50%. While there are indications fr om supernova rem- 
nant observations of one rim region (jHelder et alJl20"09l) 

2 Note that for a CR population in clusters that have been 
accelerated in structure formation shocks, the relativistic limit 
Ecn/Eth = 2(Pcii)/(.Pth) i s n ot applicable since the CR pressure 
is dominated by the transrelativistic regime. This implies a some- 
what harder equation of state for the CRs with a larger adiabatic 
index and yields the relation Bcrt/Pth = 1.6(Pgr) / {Ptb) ■ 

3 This relation should only hold for regions with long thermal 
cooling times compared to the dynamical time scale. In particular, 
it breaks down toward the center of a cooling flow cluster where 
the thermal gas cools on a shorter time scale such that the forming 
cooling flow causes adiabatic contraction of the CR population. 



as well as theoretical studies (jKang fc Jones! [2005) that 
support such high efficiencies, to date it is not clear 
whether these efficiencies apply in an average sense to 
strong collisionless shocks or whether they are realized 
for structure formation shocks at higher redshifts. Im- 
proving the sensitivity of the presented type of observa- 
tions will help in answering these profound plasma astro- 
physics questions. 

In Figure 21 we additionally compare a simulation 
where we only accelerate CRs at structure formation 
shocks with one where we additionally account for CRs 
that are injected through supernova feedback within the 
star-forming regions in our simulation. Outside the cD 
galaxy, there is no significant difference visible which sug- 
gests that the CRs injected into the ICM by supernova- 
driven winds are negligible compared with those acceler- 
ated by structure formation shocks. While this is partly 
an artifact of our simulations that neglect CR diffusion, 
we expect this behavior due to the adiabatic losses that 
CRs suffer as they expand from their compact galactic 
ISM into the dilute ICM. Assuming a conservative value 
for the density contrast of A = 10~ 3 , the CR pressure is 
diluted by P CR ~ A 4 / 3 P CR ISM ~ 10 4 PcR,ISM- 

4.3. Simplified Approach and Comparison to Previous 

Results 



As anticipated in section 13.31 there are few existing 
IACT observations of galaxy clusters; some of which 
derived limits on the CR-to-thermal pressure contained 
in clusters, in parti cular the WHIPPLE observation of 
the Perseus cluster ()Perkins et al.l I2006D and the HESS 
observations of t he A bell 85 (jAharonian" etHI I2009at 
i Domainko et al.l [2009) and Coma (jAharonian et al.l 
2009b) clusters. These works used simplifying assump- 
tions about the spectral and spatial distribution of CRs. 
They typically assumed a single CR power-law distribu- 
tion with a spectral index of T — —2.1 (that provides 
optimistic limits on the CR-to-thermal pressure) and as- 
sumed that the CR energy density is a constant fraction 
of the thermal energy density throughout the entire clus- 
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tor. Based on these two assumptions, WHIPPLE and 
HESS found in Perseus and Abell 85 E C R./E th < 0.08, 
respectively, while HESS found Ecr/E^ < 0.2 in Coma. 

To facilitate comparison with these earlier works, 
we repeated the data analysis with a spectral index 
r = —2.1 to obtain an integral upper limit .Ftjl(> 
100 GeV) = 6.22 x 10~ 12 cm~ 2 s" 1 . Following the for- 
malism of IPfrommer fc Enfilinl (|2004al) . we compute the 
7-ray flux of a CR population with T = —2.1 within 
a circular region of radius ro.6 = 0°.15 or equivalently 
200 kpc. In our isobaric model of CRs we assume that 
the CR pressure scales exactly as the thermal pressure 
and constrain EcR,/E t h < 0.053 which corresponds to 
an averaged relative pressure of {Xqr) = {Pgr) / (Pth) = 
0.033. This would be the most stringent upper limit on 
the CR energy in a galaxy cluster. 

In our adiabatic model of CRs we account for the 
centrally enhanced CR number density due to adia- 
batic contraction during the fo rmation of the cooling 
flow (jPfrommer fc Enfilinl I2004al ). We assume that the 
CRp population scaled originally as the thermal popu- 
lation but was compressed adiabatically during the for- 
mation of the cooling flow without relaxing afterward 
(we adopted tempera ture and density profiles given by 
IChurazov et al.|[2003t ). In this model, we obtain an en- 
hanced 7-ray flux level for virtually the same volume- 
averaged CR pressure or vice versa for a given flux 
limit; hence, we can put a tighter constraint on the 
averaged CR pressure. We constrain Egr/E^ < 0.03 
which corresponds to an averaged relative pressure of 
(X CR ) = (PcR>/(Pth) = 0.019. 

How can we reconcile these tighter limits with our 
simulation-based slightly weaker limit? We have to com- 
pare our simulated CR profile to a CR distribution that 
does not show any enhancement relative to the gas den- 
sity. In the central region for r < 200 kpc, we derive 
an adiabatic compression factor of 1.7 that matches that 
in our simplified approach - suggesting that our simple 
adiabatic model captures the underlying physics quite re- 
alistically. Second, we have then to relate the pressure of 
a power-law spectrum with r = 2.1 to our simulated con- 
cave spectrum. Noting that the 7-rays at 100 GeV are 
produced by CR protons at ~ 1 TeV, we normalize both 
spectra at 1 TeV and find that the simulated spectrum 
contains a larger pressure by a factor of 1.8. This factor 
brings the limit of our simplified adiabatic model into 
agreement with our simulation-based limit of the rela- 
tive CR pressure (Vcr) < 4% for the cluster core region. 
Finally, since 7-ray observations are only sensitive to the 
cluster core regions (the emission is expected to peak in 
the center due to the high target gas densities), they can- 
not constrain the average CR-to-thermal pressure within 
the entire cluster. Hence we have to use cosmological 
cluster simulations to address how much CR-to-thermal 
pressure could be additionally hidden in the peripheral 
cluster regions. 

4.4. Minimum 7-ra?/ Flux 

For clusters that host radio (mini-)halos we are able to 
derive a minimum 7-ray flux in the hadronic model of CR 
interactions. The idea is based on the fact that a steady- 
state distribution of CR electrons loses all its energy to 
synchrotron radiation for strong magnetic fields {B ^> 
Bqmb — 3.2^iG) so that the ratio of 7-ray to synchrotron 



flux becomes independ ent of the spatial distribution of 
CRs and thermal gas (|Pfrommerl [2008). This can be 
easily seen by considering the pion-decay-induced 7-ray 
luminosity L 7 and the synchrotron luminosity L v of a 
steady-state distribution of CR electrons that has been 
generated by hadronic CR interactions: 



L 1 = A 1 / dVncRn gas , 



L V = A V \ AV ucru 



e 



(a„ + l)/2 
B 



gas 



(1) 

(2) 
(3) 

Here and A v are dimensional constants that 

depend on the hadr o nic physics of the interaction 
(jPfrommer et al.l[200l iPfrommerl l200l and a v ~ 1 is 
the observed synchrotron spectral index. Hence, we can 
derive a minimum 7-ray flux in the hadronic model: 



£cmb + £b 
:A V I dVn C Rn„ as for e B > £cmb- 



T- 



7,min 



A v AitDl m 



(4) 



where L v is the observed luminosity of the radio mini- 
halo and -Dium denotes the luminosity distance to the 
respective cluster. Lowering the magnetic field would 
require an increase in the energy density of CR electrons 
to reproduce the observed synchrotron luminosity and 
thus increase the associated 7-ray flux. 

Using the values of Table [U we obtain a minimum 7- 
ray flux in the hadronic model of the radio mini-halo of 
^ r 7,min(> 100 GeV) = 6 x 10~ 13 cm -2 s -1 , assuming a 
power-law CR distribution with T ^ — 2.3. This lower 
limit is independent of the spatial distribution of CRs 
and magnetic fields. We note that the spectral index is 
consistent with the radio data 4 . It turns out that the 
requirement of strong magnetic fields violates the energy 
conditions in clusters as it implies a magnetic energy den- 
sity that is larger than the thermal energy density - in 
particular, at the peripheral cluster regions. The min- 
imum 7-ray flux condition requires a constant (large) 
magnetic field strength throughout the cluster while the 
thermal energy density is decreasing by more than a fac- 
tor of 100 from its central value. This would imply that 
the magnetic field eventually dominates the energy den- 
sity at the virial regions - a behavior that is unstable as 
it is subject to Parker-like buoyancy instabilities. Ad- 
ditionally, such a configuration would be impossible to 
achieve in first place as the magnetic energy density typ- 
ically saturates at a fixed fraction of the turbulent energy 
density which itself is only a small fraction of the thermal 
energy density in clusters (jSchuecker et al.ll2004D . Hence, 
these considerations call for lowering the assumed cluster 
magnetic fields which should strengthen the lower limits 
on the 7-ray flux considerably - however at the expense 
that these limits inherit a weak dependence on the spatial 
distribution of magnetic fields and CRs. 

Estimates of magnetic fields from Faraday rotation 
measures (RMs) have undergone a revision in the last few 

4 The CR protons responsible for the GHz radio emitting elec- 
trons are ~ 100 times less energetic than those CR protons that are 
responsible for the TeV 7-ray emission. This is consistent with the 
concave curvature found in the CR spectrum by Pinzke & Pfrom- 
mer (2010, in prep). 
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years with more recent estimates typically in the order of 
a few /iG with sli ghtly higher values up to 10 uG in cool- 
ing flow clusters (|Clarkell200 4 lEnfilin fc Vogtll2006ft . For 
the Perseus radio mini-h alo, Faraday RMs are available 
only on very small scales (jTavlor et al.l2006h , i.e. few tens 
of pc. RM estimates are of the order of ~ 7000 rad m 2 
leading to magnetic field values of ~ 25/iG assuming that 
the Faraday screen is localized in the ICM. This, how- 
ever, appears to be unlikely a s variations of 10% in the 
RM are observed on pc-scales (jTavlor et a l. 2002), while 
ICM magnetic fields are expected to be ordered on sig- 
nificantly larger scales of a few kpc (iTavlor et al.l 120061 : 
iVogt fe Enfilinl 12001 lEnfilin fe Vogtll2006D Application 
of the classical minimum-energy argument to the Perseus 
radio mini-halo data leads to estimates for the central 
magnetic field strength of B n ~ 7/xG or even Bo ~ 9/j.G 
for the more appropriate hadronic minimum-energy ar- 
gument HPfrommer fc Enfilinll2004bD . 

We select a cooling flow cluster of our sample that is 
morphologically similar to Perseus with a mass -M200 ~ 
10 15 M (D (the simulated cluster g51 of iPfrommer et ail 
2008). We adopt a conservative choice for the central 
magnetic field strength of ~ 10/xG and parameterize the 
magnetic energy density in terms of the thermal energy 
density by £b oc e^ 5 which ensures £b < £th/3 in the 
entire cluster. This allows us to strengthen the physi- 
cally motivated lower limit to J^.phys. min (> 100 GeV) = 

8.5 x 10 -13 cm -2 s _1 as shown by the dash-dotted line 
in Figure [3] In the hadronic model, this minimum 7- 
ray flux implies a minimum CR pressure relative to the 
thermal pressure. Figure [3] shows that the minimum 
flux -F-y.phys. min is a factor of 3.6 lower than the sim- 
ulated flux for Perseus in our conservative model. As 
seen in Sect. 14. 21 this model corresponds to a relative 
CR pressure of (X CR ) = (P C R.)/(Pth) = 0.04 where 
the averages represent volume averages across the entire 
cluster. Hence we obtain a minimum relative CR pres- 
sure, (A C R, mi „) = (P ca , mi„)/(Pth)/3.6 = 0.01. This 
minimum CR pressure corresponds to a minimum to- 
tal CR energy of £ C Rmin = ^CRmin/^th X E th = 

1.6 (Acr, min) x E t h = 9 x 10 61 erg where we integrated 
the te mperature and density profiles from X-ray observa- 
tions (|Churazov et al.l [2003) to obtain the total thermal 
energy of E t h = 5.7 x 10 63 erg. These considerations show 
the huge potential of combining future TeV 7-ray and ra- 
dio observations in constraining physical models of the 
non-thermal cluster emission and of obtaining important 
insights into the average distribution of cluster magnetic 
fields. 

5. DARK MATTER ANNIHILATION 

As discussed in Section [2~2l the expected DM annihi- 
lation flux is proportional to the product of a factor that 
encloses all the particle physics and a second one that 
accounts for all the involved astrophysics. Therefore, in 
order to obtain an estimate of the annihilation flux, we 
need to choose a particular particle physics model (that 
was not needed in Section 12.21 since only a comparative 
study was done there) in addition to the modeling of 
the DM distribution. Although the uncertainties in the 
particle physics factor /susy are very larg e and spread 
over s ome orders of magnitude (see, e.g., lAlbert et al.l 
2008d). it is common to use the most optimistic value 



for a given energy threshold of the telescope. This fac- 
tor just acts as a rescaling factor in the total flux, so 
we could change to the other particle physics model 
simply by rescaling for its new value. Let us assume 
/susy = 10" 32 GeV" 2 cm 3 s" 1 above 100 GeV, which 
corresponds to one of the most opti mistic allowed scenar- 
ios at the energies of interest here (jSanchez-Conde et al.1 
2007), with the neutralino as a DM particle. Then, 
taking a value of 1.4 x 10 16 GeV 2 cm -5 for the inte- 
grated astrophysical factor inside r s (as given in Sec- 
tion [22]), we obtain a maximum DM annihilation flux of 
~ 1.4 x 10 -16 cm -2 s _1 for energies above 100 GeV. The 
comparison with the derived upper limits from our obser- 
vations is not very constraining. Assuming a generic DM 
annihilation spectrum without a cutoff a nd a spectral 
index -1.5 as a good ap proximation (e.g., lAlbert et al.l 
20081 lAliuet al.ll2009ah . it can be seen from Table[2]that 
we need a boost in flux in the order of 10 4 to reach the 
predicted DM annihilation flux values, since J-ul (>100 
GeV) = 4.63 x 10~ 12 cm" 2 s" 1 . 

This boost factor could come from different uncertain- 
ties that may enhance the annihilation 7-ray flux notably 
and that were not taken into account in the above cal- 
culation. One of them, the presence of substructures, 
could play a crucial role for Perseus, as explained in sec- 
tion 12.21 Although still uncertain, its effect could en- 
hance the expected annihilation flux by m ore than a fac- 
tor of 10 for Perseus-size halos according to lKuhlen et al.l 
(2008). More recent work has shown tha t the expected 
boost factors could be as high as 200 (Springcl et al. 
I2008al| bl). However, with IACTs it is challenging to 
make use of these large boost factors as their contribu- 
tion is expected to be more important on large angular 
scales comparable to the virial extent of the cluster. De- 
tailed modeling of the substructures is needed in order 
to correctly evaluate their impact on the Perseus DM- 
induced signal. Finally, also recently proposed mecha- 
nisms in the par ticle physics side, such as the internal 
bremsstrahlung ( Brin gmann et al.l I2008D and the Som - 
merfeld effect (jLattanzi fc Silldl2009t iPinzke et al.ll2009h , 
could enhance the DM annihilation flux by more than 
one order of magnitude for some particle physics models. 

It is worth noting that the result obtained here for the 
boost factor needed in order to probe the predicted DM 
annihilation flux is comparable with prev ious observa- 
tions of the Milky Way satellite galaxies ((Albert et al.l 
I2008dt lAhu^TaTli2009al) . 

6. THE NGC 1275 EMISSION 

The SED of the NGC 1275 core is shown in Figure El 
The ra dio and optical da ta represented with gray filled 
circles (|Abdo et al.l 120091) have been obtained with low 
resolution and thus include a large contribution from the 
large-scale regions of the jet (radio) and from the host 
galaxy (optical). In the following, we model the data 
corresponding to t he core emission. This is different to 
what was done bv lAbdo et al.l (|2009l ) who used the low 
resolution data in their models. We calculated our upper 
limit, shown in Figure[5]as a red arrow, assuming that the 
spectrum in the MAGIC energy band is a power law with 
spectral index T = —2.5, as indicated by the extrapola- 
tion of the last points of the Fermi-LAT spectrum. Note, 
however, that the level of the differential upper limits is 
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only weakly dependent on the assumed spectral index 
(see Table i}. 




Log v [Hz] 



Fig. 5.— SED of the NGC 1275 core (lower two lines and data) 
and that of the well known blazar S5 0716+714 for comparison (up- 
per line and data). Gra y filled circle s are d ata points in the radio 
and optical bands from Abdo ct al. (2009). Filled black squares 
show, instead, the radio (VLBI; I Taylor et al. 1996) and the op- 
tical emission (Hubble Space Telescope (HST1; IChiaberge ct al. 
1199911 of the core alone. Th e soft X-ray bow tie is from Chan- 
dra iBalmavcrdc ct al. 2000!) , while t he red filled circles represent 
the Fermi-LAT spectrum taken from Abdo ct al. (2009]). The red 
arrow shows the MAGIC upper limit between 80 and 100 GeV. 
The lower blue and red lines show the emission from the spine and 
the layer of the structured jet, respectively. The upper blue lines is 
the SED of the spine as observed at a small angle (see the text for 
details); fo r comparison, we report hist orical data of S5 0716+714 
(data from Tavccchio & Ghiscllini 2009 and references therein). 

The data clearly show a double-peak SED. The radio- 
optical data suggest a peak of the emission in the IR 
band, simila r to the case of other 7-ray emit ting ra- 
dio galaxies (|Tavecchio fc Ghisellinil [20081 [20091) . High- 
energy data constrain the peak frequency of the sec- 
ond component at about 1 00 MeV. As discussed in 
iTavecchio fc Ghisellinil (|2008t ). in these cases a one-zone 
synchrotron self-Compton (SSC) model for the entire 
emission is implausible, since the large separation in fre- 
quency between the two peaks would require extremely 
large values of the Doppler factor: 

S ~ 258 L sA2 .sL c ^l 4 v~i 3 Q ^c,23 Rm ( 5 ) 

where L s — v a L(y s ), Lq — v c L(v c ), v s and vq are the 
synchrotron and SSC peak luminosities and frequencies, 
respectively, and R is the size of the emitting region. 
Here Q = W X Q X in cgs units, and we use the values 
derived for NGC 1275. In this estimate, we assume the 
typical size of the emission regions derived in blazars, 
R = 10 16 cm, though the Fermi-LAT data do not al- 
low us to constrain the radius of the e mission region us- 
ing the variability (|Abdo et al.l 120091 ) . Such large val- 



ues of the Doppler factor are rather unlikely. Typical 
values found modeling the SED of blaz ars are around 
10 - 20 (e.g., iCelotti fc Ghisellinil [2001 . with few ex- 
treme TeV BL Lacs requiring larg er values during excep- 
tional states (5 ~ 50 - 100; e.g., IBegelman et al.l 120081: 
IGhisellini fc Tavecchiol [20081 ) . Arguments based on the 
observation of superluminal motions a t Very Long Base- 
line I nterferometry (VLBI) scales (e.g.. lKellermann et al.l 
120041 ) and on the unification of blazars with radio galax- 
ies al so require values around 10 (e.g. JUrrv fc Padovanil 
H995ft . 

The most direct way to overcome the problem posed 
by the large Doppler factor is to assu me two emission 
region s, as in the spine-layer model of Ghi sellini et al.l 
(2005). In this scenario the jet is assumed to be struc- 
tured, with a fast inner region (the spine) surrounded 
by a slower sheet (the layer). Both components produce 
synchrotron and IC radiation and they are in radiative 
interplay: the synchrotron radiation from one component 
is seen boosted (by the relative velocity) by the other one 
and thus the IC emission of both regions is enhanced with 
respect to the standard SSC. In radio galaxies, in which 
the jet is observed at relatively large angles, the emission 
is expected to be dominated by the layer, which, due to 
the lower bulk Lorcntz factor, has a larger emission cone. 
At a smaller angle, instead, the emission is dominated by 
the spine, as in blazars. 

We reproduce the SED with the spine-layer model. 
The orange line in Figure 03 shows the emission from the 
layer, while the spine produces the emission shown by the 
blue bottom line. The spine is assumed to be a cylinder 
of radius R = 1.5 X 10 16 cm, height Hs = 1.5 x 10 16 cm 
(as measured in the spine frame), and in motion with 
bulk Lorentz factor Tg = 15. The layer is modeled as an 
hollow cylinder with internal radius R, external radius 
i?2 = 1.2 x R, height = 4 x 10 16 cm (as measured in 
the frame of the layer), and bulk Lorentz factor Ll = 3. 
Each region contains a tangled magnetic field with inten- 
sity B§ = 2.5 G and £?l = 1 G, and it is filled by relativis- 
tic electrons assumed to follow a (purely phenomenolog- 
ical) smoothed broken power-law distribution extending 
from 7 m j n to 7 max and with indices m, n 2 below and 
above the break at 7b. For the spine we use 7 m i n = 40, 
7b = 2 x 10 4 , 7 max = 10 5 , n\ — 2, n 2 = 3.5. For the layer 
7min = 10, 7b = 4 x 10 3 , 7 max = 10 5 , m = 2.4, n 2 = 4.2. 
The normalization of these distributions is calculated as- 
suming that the systems produce an assumed (bolomet- 
ric) synchrotron luminosity L' s = 10 42 erg s _1 and 
L' syn L = 2.7 x 10 41 erg s _1 (as measured in the local co- 
moving frame of the spine and layer, respectively), which 
is an input parameter of the model. As said above, the 
seed photons for the IC scattering are not only those pro- 
duced locally in the spine (layer), but we also consider 
the photons produced in the layer (spine). We assume a 
viewing angle of 9 = 15°. As discussed above, the same 
jet observed at a smaller angle would be dominated by 
the emission from the spine and we expect that its SED 
resembles those of typical blazars. We show the SED of 
the jet when observed at an angle of 4°. 5 (blue upper 
line in Figure 03) . The SED is dominated by the emis- 
sion of the spine. For comparison, we report historical 
data for the well-known blaz ar S5 0716+714 (data from 
ITavecchio fc G hiscllini 2009] and references therein) . 
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Note that, as observed, the model naturally predicts a 
very rapid decrease of the emission level above 10 GeV, 
due to the decreasing efficiency of the IC scattering in 
the Klein-Nishina regime. The position of this break 
is tightly constrained by the Fermi-LAT spectrum and 
MAGIC upper limit. In our model, this is critically de- 
pendent on the value of the frequency of the target pho- 
tons for the IC scattering that in the spine-layer scenario 
are mainly those coming from the spine (and scattered 
by the electrons in the layer). Therefore, the determina- 
tion of the cut-off frequency between the F ermi-L AT and 
the MAGIC band allows us to infer the peak frequency 
of the synchrotron component of the spine. For instance, 
assuming that the peak of the spine is at IR frequencies or 
below (using for the layer the same parameters adopted 
above), we predict a flux in the MAGIC band above the 
measured upper limit. This argument allows us to fix 
the synchrotron peak of the spine at optical- UV frequen- 
cies. This, in turn, assures that the beamed counterpart 
of NGC 1275 is an intermediate BL Lac object, as the 
chosen S5 0716+714. In conclusion, the knowledge of the 
upper limit at the low-energy end of the MAGIC band 
offers us the important possibility of having independent 
limits on the characteristics of the emission of the (oth- 
erwise invisible) spine and thus of constraining the kind 
of beamed counterpart of this radio galaxy. Future ob- 
servations can confirm or rule out our interpretation. In 
particular, the detection of photons above <~ 100 GeV 
would be challenging for the scenario depicted here, re- 
quiring major changes in the emission properties of the 
spine. 

7. CONCLUSIONS 

The Perseus cluster was observed by MAGIC during 
2008 November and December resulting in 24.4 hr effec- 
tive observation time of very high data quality. No sig- 
nificant excess of 7-ray was detected beyond the energy 
threshold of 80 GeV. 

Using simplified assumptions (power-law CR spectra, 
constant ratio of CR-to-thermal energy density) that 
have been adopted in earlier work, we obtain a limit 
on the CR energy of Ecu/E t h < 5%. This limit could 
be tightened furthermore by considering an adiabatically 
contracted CR population during the formation of the 
cooling flow yielding EcR/E t h < 3%. This would be the 
most stringent constraint on the CR energy using 7-ray 
observations to date. Using cosmological cluster simula- 
tions, it turns out that these assumptions are not fulfilled 
for CR populations that have been accelerated by struc- 
ture formation shocks: while the adiabatic model seems 
to match the simulated CR profiles toward the center 
very well, the expected ratio of CR-to-thermal pressure 
increases toward the peripheral cluster regions causing 
the volume-averaged pressure across the entire cluster to 
increase by a factor of 2. In addition, the CR spectral 
distribution shows a concave curvature with a spectrum 
that flattens toward high energies with a spectral index 
of T ~ -2.2 in the TeV regime. This implies that the CR 
pressure is enhanced by an additional factor of almost 2. 
Using our simulated flux, we obtained an upper limit on 
the CR-to-thermal pressure averaged across the entire 
cluster volume of (Xqk) < 8% and < 4% for the cluster 
core region. This corresponds to an upper limit on the 
CR energy of Ecr/E^ < 13% and < 6.5%, respectively. 



We note that this is the first work where results from 
cosmological simulations and observational data analysis 
are combined. This demonstrates the need for cosmo- 
logical simulations in order to more reliably predict CR 
spectra which provide a safeguard against too simplified 
and optimistic models which then lead to limits that are 
too tight. 

The upper limits resulting from the data analysis are 
a factor of ~ 2 larger than our conservative model pre- 
diction for the CR-induced 7-ray emission and hence in 
agreement with our cosmological cluster simulations. Fu- 
ture more sensitive measurements will be able to put 
interesting constraints on the maximum shock accelera- 
tion efficiency. Using minimum 7-ray flux arguments, we 
show that improving the sensitivity of this observation by 
a factor of about 7 will enable us to finally critically test 
the hadronic model for the Perseus radio mini-halo: a 
non-detection of 7-ray emission at this level implies CR 
fluxes that are too small to produce enough electrons 
through hadronic interactions with the ambient gas pro- 
tons to explain the observed synchrotron emission. 

As DM dominates the cluster mass, significant 7-ray 
emission resulting from its annihilation is also expected. 
With the assumed particle p hysics model, one of the most 
optimistic allowed scenarios ([Sanchez-Conde et al.|[2007n 
with the neutralino as a DM particle, the boost factor 
for the typically expected DM annihilation-induced emis- 
sion is constrained to < 10 4 . Note that for this estima- 
tion, we neglected possible contributions from internal 
bremsstrahlung, Sommerfeld enhancement as well as en- 
hancement factors due to substructures. 

The upper limits obtained for the NGC 1275 emission 
are consistent with the recent detection by the Fermi- 
LAT satellite. In this case a one-zone SSC model for 
the entire emission is implausible, since the large sepa- 
ration in frequency between the two peaks would require 
extremely large values of the Doppler factor for the jet 
(|Tavecchio fc Ghiseliim] 120081 ) . The most direct way to 
overcome this problem is to assume two emission r egions, 
as in the spine-layer model (iGhisellini et a l. 2005) which 
explains the radio galaxy emission. 

While no galaxy cluster has been detected in 7-rays up 
to now, our estimations indicate that Perseus is among 
the most promising clusters to be detected by IACTs. 
Using the newly inaugurated MAGIC seco nd telescope 
and o perating the telescopes in stereo mode (|Colin et al.l 
2009), a total observation time of about 150 hr may give 
us a chance to detect the CR-induced 7-ray emission or 
to definitively probe the validity of the hadronic model of 
radio (mini-)halos. As the emission of NGC 1275 domi- 
nates the accessible energy range of the Fermi-L AT satel- 
lite, it could potentially hinder the satellite from detect- 
ing the CR as well as the DM-induced 7-ray emission in 
this cluster. Similar problems might arise in other clus- 
ters. Therefore, the IACTs will play a crucial role in the 
quest for 7-ray emission from galaxy clusters. 
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